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A beam of polar molecules can be focused and transported through an ac electric quadrupole guide. 
At a given ac frequency, the transmission of the guide depends on the mass-to-dipole-moment (m/p) 
ratio of the molecular quantum state. Here we present a detailed characterization of the m/ji selector, 
using a pulsed beam of benzonitrile (CeHsCN) molecules in combination with rotational quantum 
state resolved detection. The arrival time distribution as well as the transverse velocity distribution 
of the molecules exiting the selector are measured as a function of ac frequency. The pi/ Api resolution 
of the selector can be controlled by the applied ac waveforms and a value of up to 20 can be obtained 
with the present setup. This is sufficient to exclusively transmit molecules in the absolute ground 
state of benzonitrile, or rather in quantum states that have the same m/p value as the ground state. 
The operation characteristics of the m/ji selector are in quantitative agreement with the outcome 
of trajectory simulations. 
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I. INTRODUCTION 

Molecular beams have played central roles in many ex- 
periments in physics and chemistry and have found a 
wide range of applications [1 . For many of these appli- 
cations it is crucial that the internal quantum state distri- 
bution of the molecules in the beam is accurately known 
and can be controlled. Cooling of the molecules in an 
adiabatic expansion strongly reduces the number of pop- 
ulated levels, thereby tremendously simplifying the opti- 
cal spectra of the molecules, for instance [7. For large 
molecules, however, many levels will still be populated 
at the low temperatures (~1 K) that can be achieved in 
a molecular beam. Moreover, even at the lowest possi- 
ble temperatures, molecules with a complex potential en- 
ergy landscape like bio- molecules, are known to be frozen 
into different minima and to be present in various con- 
formational structures in the beam [3]. To control the 
quantum-state distribution (and thereby the conforma- 
tional distribution) of large molecules in a beam, there- 
fore, additional filtering techniques need to be applied. 
This can be achieved with electric, magnetic or optical 
fields, or any combination thereof [4]. In the work pre- 
sented here we discuss the quantum state selection of a 
beam of large molecules using electric fields. 

A wide variety of electric field geometries has been 
used in the past to manipulate the trajectories of polar 
molecules in a beam. The original geometries were de- 
vised to create strong electric field gradients on the beam 
axis to efficiently deflect molecules. Later, electric field 
geometries were designed to focus molecules in selected 
quantum states further downstream. For these experi- 
ments, one has to distinguish between molecules that are 
in so-called low-field-seeking or high-field-seeking quan- 
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tum states; as the name suggests, the former are at- 
tracted to regions of low electric field whereas the latter 
experience a force in the direction of a high electric field. 
In a seminal experiment, performed more than half a cen- 
tury ago, an electrostatic quadrupole focuser was used to 
couple a beam of ammonia molecules into a microwave 
cavity. This beam initially contained an almost equal 
number of molecules in the upper (low-field-seeking) and 
the lower (high- field-seeking) level of the | J^K) =| 1, 1) 
inversion doublet. The radially symmetric electric field 
in the quadrupole focuser had a zero of the electric field 
strength on the beam axis, resulting in a focusing (defo- 
cusing) of ammonia molecules in the upper (lower) inver- 
sion doublet level. The resulting inverted population dis- 
tribution in the microwave cavity led to the demonstra- 
tion of microwave amplification by stimulated emission 
of radiation, i.e., the invention of the maser, by Gordon, 
Zeiger and Townes |5]. 

Focusing is considerably more complicated for 
molecules in high-field-seeking states than for molecules 
in low-field-seeking states as no electric field maximum 
can be created on the molecular beam axis using static 
fields alone. Due to their small rotational constants 
and the resulting high density of levels, however, large 
molecules are high-field seeking in all their low-lying ro- 
tational levels already at modest electric field strengths. 
To focus large molecules, therefore, dynamic focusing has 
to be applied [6 , much like what is routinely done to 
transport charged particles . Dynamic focusing of neu- 
tral polar molecules has first been experimentally demon- 
strated on ammonia molecules in high-field-seeking quan- 
tum states [8 . More recently, it has been applied in 
the deceleration of metastable CO [9] in the trapping of 
ND3 molecules [10 , in the guiding of a cold beam of 
molecules from an effusive source [11 and in the trans- 
port of CaF [12], for instance. Dynamic focusing has also 
been applied to large molecules, e.g., in the deceleration 
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of benzonitrile fT3] as well as in the separation of two 
different conformational structures of aminophenol p!4] . 
In the latter experiment the ac electric quadrupole guide 
was used to selectively transmit molecules depending on 
their mass-to-dipole-moment (m/p) ratio, analogous to 
the transmission of ions with a certain mass-to-charge 
(m/q) ratio in a quadrupole mass filter. 

Here we present a detailed characterization of a second 
generation m/p selector. A pulsed beam of rotationally 
cold benzonitrile (CeHsCN) molecules is injected into the 
ac electric quadrupole selector, and the arrival time dis- 
tribution as well as the transverse velocity distribution 
of the molecules exiting the selector are measured using 
high-resolution electronic spectroscopy. The dependence 
of the p/Ap resolution of the selector on the applied ac 
waveforms is studied, and the transmission of selected 
subsets of rotational levels is demonstrated. The opera- 
tion characteristics of the m/p selector are compared to 
the outcome of numerical trajectory simulations. 



II. EXPERIMENTAL SETUP 



(a) 



A scheme of the experimental setup is shown in |Fig-| 
ure l[ a). Argon gas at 2.5 bar is bubbled through liq- 
uid benzonitrile in a room-temperature container. The 
resulting gas mixture of 0.04% benzonitrile in Argon is 
expanded through a 0.8 mm diameter nozzle into vacuum 
using a pulsed valve (General Valve, Series 99) at a rep- 
etition rate of 40 Hz. A beam of benzonitrile molecules 
with a rotational temperature of about 1 K and with a 
mean forward velocity of around 570 m/s is thus pro- 
duced. The value for the most probably velocity is not 
precisely known but is expected to be close to, but prob- 
ably slightly higher than, that for a mono-atomic ideal 
gas expanding from a room-temperature valve [1 ; with 
Vp = ^/5kBTo/m with ks the Boltzmann constant, m the 
mass of the Ar atoms and To=295 K, a lower limit of Vp 
= 552 m/s is obtained. For a detailed analysis of the ex- 
perimental results, the initial length of the gas pulse and 
the width of the velocity distribution are other impor- 
tant parameters. We do not know these independently, 
but we can accurately determine the relative spread in 
the arrival time distribution at the end of the beam ma- 
chine as being 4.8% {vide infra). We mainly attribute 
this to the width of the velocity distribution, which is as- 
sumed to have a full-width-at-half maximum (FWHM) of 
about 25 m/s, corresponding to a relative velocity spread 
of about 4.4%. This narrow velocity distribution cor- 
responds to a longitudinal translational temperature of 
0.5-0.6 K. 

After passing through a 1.5 mm diameter skimmer 
(Beam Dynamics, model 2) placed about 40 mm from 
the nozzle, the molecules enter the selector 32 mm be- 
hind the tip of the skimmer. The selector consists of four 
1.64 m long and 6 mm diameter stainless steel electrodes. 
The ends of the highly polished electrodes are rounded 
off with a radius of curvature of 3 mm. The electrodes 
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FIG. 1. (a) Schematic view of the experimental setup, (b) 
Photograph of the end of the selector, showing how the cylin- 
drical electrodes are mounted to the four rigid backbones, 
(c) Electric field for either one of the two high-voltage con- 
figurations. The electric field at the center is 48 kV/cm and 
contour lines are drawn every 10 kV/cm away from the center. 
The transparent cross-shaped overlay shows the area in the 
(x,y)-plane to which the molecules are confined, (d) Switching 
scheme for two different duty cycles at the same ac frequency 
(l/T). 



are centered at the corners of a square with a diagonal 
of 10.5 mm. This implies that the closest distance of the 
surface of each of the four electrodes to the molecular 
beam axis is 2.25 mm; the closest distance b etween ad ja- 
cent electrodes is 1.42 mm, as indicated in Figure l| ^b). 
To achieve the best possible performance of the selector 
care is taken to keep the distances between the electrodes 
as constant as possible over the full length of the device. 
To this end, four rigid 1.4 m long aluminum holders have 
been machined that serve as straight, stable backbones 
on which the electrodes are mounted. The electrodes, 
which are out of magnetic stainless steel, are attached 
via small magnets that are inserted from the backside 
into the holders. This mounting scheme prevents the 
unavoidable bending of the electrodes that would occur 
when screw holes would be drilled into them and it also 
prevents mechanical stress that would be present when 
the electrodes would be screwed onto the holders. In the 
completely mounted setup, we can only accurately mea- 
sure the spacing between adjacent electrodes. It is found 
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that this spacing deviates less than ±0.02 mm from its 
nominal value of 1.42 mm over the full length of the se- 
lector. 

High voltages are connected to the four electrodes of 
the selector according to either one of the configurations 
shown in Figure ij ^c). In configuration 1, voltages of 
+9.5 kV are applied to the upper two electrodes while 
the lower two electrodes are at -9.5 kV, creating a saddle- 
point of the electric field of 48 kV/cm on the molecu- 
lar beam axis. In the horizontal plane, the electric field 
increases with distance from the center whereas it de- 
creases away from the center in the vertical plane. Ben- 
zonitrile molecules in high-field-seeking states, therefore, 
will be focused in the vertical direction while being de- 
focused h orizontally, as indicated by the white arrows in 
Figure 1 'c) . In going from configuration 1 to configura- 



tion 2, the voltage applied to the right upper (left lower) 
electrode is switched from +9.5 kV to -9.5 kV (-9.5 kV 
to +9.5 kV), effectively rotating the field configuration, 
and thereby the (de-)focusing direction, by 90 degrees. 
By switching back and forth between these two config- 
urations, benzonitrile molecules can be transported on 
stable trajectories through the selector. The switching 
between the configurations is performed with two push- 
pull high-voltage switches (Behlke, HTS 301- 03-GSM) 
according to the scheme depicted in [Figure l| d) . The 
time-dependence is shown there as a block-function; in 
reality it takes about 0.5 ps to reach the final voltage. 
The ac switching frequency (1/T) is in the 2-8 kHz range 
and the fraction of the period during which configuration 
1 is applied is denoted as the duty cycle d. By vary- 
ing the duty cycle the p/Ap resolution is changed [T5 
and in the present experiments values of d in the range 
0.5-0.6 have been used. The cross-shaped overlay in 



Figure 1 ^c) depicts the typical area in the (x,y)-plane 
to which molecules stay confined while fiying in the z- 
direction. From this, benzonitrile molecules transmit- 
ted on stable trajectories are seen to sample electric field 
strengths in the range of 25-77 kV/cm. 

The benzonitrile molecules are monitored 40 mm be- 
hind the end of the selector via laser induced fiuorescence 
(LIF) detection. Individual rotational transitions in the 
origin band of the electronic ^ Sq system are induced 
using narrowband radiation from a continuous wave ring 
dye laser (Coherent 899-21) that is frequency doubled in 
an external cavity (Spectra Physics, LAS Wave Train). 
Typically 15 mW of tunable 274 nm radiation with an 
instantaneous bandwidth of less than 1 MHz is obtained 
in this way. The laser frequency is actively stabilized 
against a frequency stabilized Helium-Neon laser (SIOS 
SL 03), thus allowing to keep the frequency constant to 
within 2 MHz for several hours. The collimated 2.5 mm 
diameter laser beam intersects the molecular beam under 
a right angle in the horizontal plane while the fiuores- 
cence is detected from above, i.e., along the y-direction. 
The transmission of benzonitrile molecules through the 
selector in individual rotational JxaK^ levels is measured 
via time-resolved counting of the LIF photons. In this 
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FIG. 2. Rotationally resolved fluorescence excitation spectra 
of benzonitrile. The horizontal axis is relative to the fre- 
quency of the 111 ^ Ooo transition at 36512.74 cm~^. Middle 
trace: Spectrum obtained without voltages applied to the se- 
lector, i.e., in free flight. Lower trace: Simulated free- flight 
spectrum (plotted upside down), assuming a rotational tem- 
perature of 1.2 K. The intensities are scaled to the intensity of 
the 111 ^ Ooo transition in the middle trace and a Voigt pro- 
file with a Lorentzian contribution (FWHM) of 8 MHz and a 
Gaussian contribution (FWHM) of 3.7 MHz is taken. Upper 
trace: Spectrum obtained with the selector operating at an ac 
switching frequency of 3.8 kHz, (i=0.50. The J^/ j^, ^ ^KaKc 
assignment for each of the rotational lines is indicated. 



setup, we are also sensitive to the transverse velocity dis- 
tribution of the benzonitrile molecules in the x-direction 
via the Doppler-profile of the rotational transitions. 



III. EXPERIMENTAL RESULTS 



In Figure 2 a selected part of the rotationally resolved 
fluorescence excitation spectrum of benzonitrile is shown. 
The fluorescence signal is integrated over a time-interval 
of 0.5 ms, thereby accumulating the signal from all ben- 
zonitrile molecules in a pulse. The spectrum shown in 
the middle trace is recorded in free flight, i.e., with- 
out voltages applied to the selector. A total of nine 
clearly separated rotational lines, originating from nine 
distinct low-lying rotational levels, are recognized. From 
this spectrum, the rotational temperature of benzonitrile 
in the molecular beam can be accurately determined. 
Taking the known molecular constants of benzonitrile in 
the ground state [16 and in the electronically excited 
state [17 in combination with the degeneracy of the lev- 
els (spin statistics) and the rotational line intensities in 
the electronic transition , the observed spectrum can 
be simulated with the rotational temperature as the only 
adjustable parameter. A simulated spectrum for a rota- 
tional temperature of 1.2 K is shown in the lower trace, 
and is seen to match the experimental spectrum rather 
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3. Expanded view of the In ^ Ooo transition shown 



Figure 2 measured with the selector on (upper trace) and 
off (lower trace). Voigt profiles fitted to the experimental 
data are shown and the widths (FWHM) of the Lorentzian 
and Gaussian contributions to the Voigt profiles are given. 
The difference of the observed and fitted profiles is shown 
underneath, slightly offset for clarity. 



well. A Voigt profile with a Lorentzian (Gaussian) con- 
tribution with a width (FWHM) of 8 MHz (3.7 MHz) 
is taken in the simulation. The width of the Lorentzian 
contribution has been taken from earlier high-resolution 
spectroscopy experiments [17 and is consistent with our 
present observations. It should be noted that the radia- 
tive lifetime of the vibrationless level in the Si state is 
about 70 ns [I9l[20|, thus contributing less than 2.5 MHz 
to the observed homogeneous linewidth; the remaining 
contribution has to result from pure dephasing, with an 
associated pure dephasing time of about 60 ns. In the up- 
per trace, the spectrum obtained with the ac quadrupole 
selector operating at a switching frequency of 3.8 kHz, 
and with a duty cycle of 0.50, is shown. The fluores- 
cence intensity is plotted on the same scale as for the 
spectrum recorded in free flight and an increase in peak 
intensity of almost a factor five is obtained for benzoni- 
trile molecules in their rotational ground state (Ooo) level. 
For rotational transitions starting from rotationally ex- 
cited levels the intensity increase upon using the selector 
is seen to be smaller, however, and the intensity even de- 
creases for the transitions starting from the 5o5 and the 
6o6 levels. 

It is evident from the spectra shown in [Figure 2| that 
the rotational transitions are broadened when the selec- 



tor is switched on. This is seen more clearly in Figure 3 
where an expanded view of the In ^ Ooo transition, 
recorded with the selector on (upper trace) and off (lower 
trace), is shown together with the best fitting Voigt pro- 
files. As the Lorentzian contribution to the profiles is 
kept fixed to 8 MHz, the width of the Gaussian contri- 
bution and the peak intensity are the only fitting param- 
eters. A Gaussian contribution with a width of 3.7 MHz 



is obtained under free-fiight conditions, determined by 
the residual Doppler broadening and the bandwidth of 
the laser system. When the selector is on, the width 
of the Gaussian contribution increases to 19.1 MHz due 
to the increased transverse velocity of the benzonitrile 
molecules exiting the selector. The spectrally integrated 
line intensity increases by a factor of 9.3 when the selector 
is switched on. As the transverse spatial distribution of 
the molecules exiting the selector will be larger, the over- 
lap with the detection laser will be somewhat smaller, 
implying that this factor 9.3 is a lower limit for the over- 
all increase of the transmission of benzonitrile molecules 
in the Ooo level upon switching on the selector. 
The 
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spectral line-profiles shown in 
recorded by accumulating the fluorescence signal from all 
ground-state benzonitrile molecules in the pulse. When 
these line profiles are detected for narrow time intervals, 
the time dependence of the transverse velocity distribu- 
tion of the molecules exiting the selector can be mea- 
sured. Such measurements are shown in two-dimensional 



contour-plots in Figure 4 For this, the intensity of the 
fluorescence signal, binned in 2 ps time-intervals, is plot- 
ted in a false-colour representation as a function of the 
frequency of the excitation laser, which is scanned in 
1 MHz steps from -40 MHz to +40 MHz around the cen- 
ter frequency of the In ^ Ooo transition. Interesting 
structures are observed in these two-dimensional plots, 
that are clearly different when the selector is operated at 
3.45 kHz (left contour plot) or at 3.55 kHz (right contour 
plot). The panel on the left shows the time-integrated 
LIF signals as a function of laser frequency. The Doppler 
shift on the vertical axis can be directly translated into 
a velocity component along the propagation direction of 
the laser, i.e., in the x-direction, and the corresponding 
velocity scale is given as well. The time on the horizontal 
axis is the time (in ms) since the triggering of the pulsed 
valve; the actual opening of the pulsed valve has a delay 
with respect to the trigger of approximately 0.2 ms. The 
panels at the bottom show the arrival time distributions 
integrated over the Doppler-broadened profiles that are 
seen to have a FWHM width of about 0.15 ms. 

In the experiments the switching sequence is started 
prior to the molecules entering the selector. The trans- 
verse velocity distribution that the molecules have in the 
detection region is mainly determined by the forces that 
the molecules experience near the end of the selector. For 
this reason, the phase in the switching cycle at which 
the molecules exit the selector, referred to hereafter as 
the end phase, is schematically shown above the contour- 
plots. As the detection of the molecules is performed only 
40 mm from the end of the selector and as the width of 
the longitudinal velocity distribution is sufficiently nar- 
row, the arrival time of the molecules in the detection 
region can be directly mapped to an end phase. 

In the measurements recorded with an ac switching 
frequency of 3.45 kHz, the molecules arriving at 3.3 ms 
exited the selector in the middle of configuration 2, i.e., 
during the middle of a horizontal focusing phase. This 



5 



^ 




3.2 3.4 
Intens. (arb.u.) Arrival time (ms) 

1.0 



3.2 3.4 
Arrival time (ms) 




FIG. 4. Intensity of the LIF signal, plotted in a false-colour 
representation, as a function of the frequency of the excitation 
laser and as a function of arrival time for two different ac 
switching frequencies. The vertical axis is also given in terms 
of the velocity in the x-direction. Above the contour-plots, the 
phase in the switching cycle at which the molecules exited the 
selector is indicated. The panel on the left shows the time- 
integrated LIF signals as a function of Vx- The panels at the 
bottom show the arrival time distributions integrated over the 
Doppler-broadened profiles. 
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is actually at a turning point of the micro motion in the 
selector, and thus results in a rather parallel beam and a 
correspondingly narrow spectral line-profile. In the mea- 
surements recorded at the slightly higher frequency of 
3.55 kHz, the molecules arriving at 3.3 ms exited the se- 
lector at the end of a horizontal focusing phase. These 
molecules therefore have the maximal velocity compo- 
nents along the x-direction, leading to the broadest line- 
profile; benzonitrile molecules exiting the selector at the 
middle of either configuration 2 or configuration 1 again 
have a narrow spectral profile. Related experiments to 
visualize the motion of dynamically confined neutral par- 
ticles have been performed on metastable CO molecules 
passing through an alternating gradient decelerator ^ 
as well as on Rb atoms in an ac electric trap [22] . 



In the upper panel of Figure 5 , the time-integrated 
LIF signal of benzonitrile molecules in the Oqo level is 
shown as a function of the ac switching frequency in the 
2-8 kHz range. The ac frequency is scanned in steps of 
50 Hz for all curves except for the one at a duty cycle of 
0.59, which is recorded in 10 Hz steps; the latter curve 
is shown once more on an expanded scale in the inset. 
The excitation laser is kept fixed to the center of the 
111 ^ Ooo transition and the signal intensity is normal- 
ized to the signal that is obtained when no high voltages 
are applied to the selector. To have well-defined condi- 
tions in the experiment, the timing of the ac switching 
sequence was arranged such that 0.2 ms after the trigger 



FIG. 5. Upper panel: Time-integrated LIF signal intensity as 
a function of the ac switching frequency. Transmission curves 
are measured at six different duty cycles, indicated in the in- 
set, while the laser frequency is kept fixed to the center of the 
111 ^ Ooo transition. The transmission curve for (i=0.59 is 
also shown on an expanded scale. The signal is normalized to 
the LIF intensity obtained in free flight. Lower panel: Mea- 
sured intensity of the LIF signal, plotted in a false-colour rep- 
resentation, as a function of the ac switching frequency and as 
a function of arrival time for three different duty cycles. The 
panel on the left shows the corresponding time-integrated LIF 
signals as a function of ac switching frequency. 



to the pulsed valve, i.e., at the moment that the pulsed 
valve opens, there is a fixed phase for all ac frequencies. 
Although this implies that the phase that a benzonitrile 
molecule experiences at the entrance of the selector, the 
so-called start phase, is different for different ac frequen- 
cies this start phase can be precisely calculated once the 
position of the molecule in the pulse and its velocity is 
known. The end phase is in this situation different for 
different ac switching frequencies as well, leading to the 
observed intensity modulation on the individual curves. 
All the different curves show no transmission of the selec- 
tor below a certain cut-off frequency, followed by a steep 
onset of the transmission around 3 kHz, a maximum at 
frequencies between 3-4 kHz and a gradual decrease in 
transmission for higher frequencies. This general behav- 
ior has been described before [TO^, and is similar to a 
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typical transmission curve in a quadrupole mass filter, for 
instance. At low ac switching frequencies the trajectories 
of the molecules are unstable, whereas at high frequen- 
cies the effective potential well in which the molecules 
are confined becomes very shallow. Only in an interme- 
diate frequency range good transmission on stable trajec- 
tories is achieved. The transmission of the m/p selector 
can be strongly modified when an asymmetric switch- 
ing cycle is applied, i.e., when values of d ^ 0.50 are 
chosen [15]. Transmission curves for six different duty 
cycles are shown to illustrate this. For a duty cycle of 
0.50 the FWHM of the transmission curve for benzoni- 
trile molecules in their ground-state level is about 1 kHz. 
With increasing duty cycle, the width of the transmis- 
sion curve gets considerably narrower while the trans- 
mission intensity goes down. At a duty cycle of d = 0.59, 
the FWHM of the transmission curve has come down to 
0.16 kHz, while the peak transmission of the selector is 
then reduced by about a factor ten. 



In the lower panel of Figure 5 , the intensity of the mea- 
sured LIF signal is shown in a false-colour representation 
as a function of ac switching frequency and as a func- 
tion of arrival time for three selected duty cycles. The 
transmission curves shown in the upper panel have been 
obtained by time-integrating these signals, and three of 
these are shown again in the lower panel on the left. To 
understand the two-dimensional plots, one has to realize 
that at a given ac frequency and with the excitation laser 
at the center of the rotational transition, thereby only 
probing molecules with small velocities in the x-direction, 
the most strongly peaked arrival time distribution is ob- 
tained for molecules exiting the selector in the middle of 
configuration 2. This is directly evident from inspecting 
the contour-plots shown in [Figure^ As shown there, at 
an ac frequency of 3.45 kHz, these benzonitrile molecules 
produce a narrow arrival time distribution, peaked at 
3.3 ms. With increasing ac frequency, the correspond- 
ing peak in the arrival time distribution shifts to earlier 
times, explaining the slopes of the intensity profiles ob- 
served in [Figured When the molecules exit the selector 
during the middle of configuration 1, i.e., during a verti- 
cal focusing phase, the overall intensity of the LIF signal 
is less. This can be understood from simple arguments: 
at the middle of either focusing phase, the beam coming 
out of the selector will be most parallel and its spatial 
distribution at the end of the selector will be maintained 
in the detection region. In the middle of a vertical fo- 
cusing phase (configuration 1) the cloud of molecules is 
necessarily larger in the y-direction than in the middle 
of a vertical defocusing phase (configuration 2); this is 
actually the underlying operation principle of alternat- 
ing gradient focusing. The spatial overlap with the laser 
beam will thus be smaller in the middle of configuration 
1. Moreover, the velocity spread in the horizontal di- 
rection is larger in the middle of configuration 1 than in 
the middle of configuration 2. Both the decreased over- 
lap in the vertical direction and the increased velocity 
spread in the horizontal one lead to a lower LIF intensity 
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FIG. 6. Rotationally resolved fluorescence excitation spec- 
tra of be nzonitrile. Upper trace: Free- flight spectrum from 
Figure 2 in which the J^/ ^, ^ JKaKc assignment for each 



or the rotational lines is indicated. Middle trace: Spectrum 
obtained with the selector operating at an ac switching fre- 
quency of 2.7 kHz, (i=0.50. Lower trace: Spectrum obtained 
with the selector operating at an ac switching frequency of 
3.8 kHz, (i=0.565. 



in the middle of configuration 1 than in the middle of 
configuration 2. This explains the occurrence of maxima 
with alternating intensity with increasing ac frequency, 
for a given arrival time. Of course, the profile of the ar- 
rival time distribution at any ac frequency is weighed by 
the intensity profile of the pulse, which is peaked around 
3.3 ms and has a (FWHM) width of 0.15 ms. It is clear 
from these two-dimensional plots that when these are in- 
tegrated over time to obtain the transmission curves as 
function of the ac frequency, some remnants of the struc- 
ture originating from the varying end phases remain. 



The transmission curves shown in Figure 5 are for ben- 
zonitrile molecules in the Oqo ground-state level. The 
overall transmission profile shifts to lower ac switching 
frequencies when the forces on the molecules are weaker, 
i.e., for benzonitrile molecules in rotational levels with a 
lower dipole moment. This feature can be used to selec- 
tively transmit benzonitrile molecules in s pecific (su bsets 
of) rotational quantum levels, as shown in Figure 6 The 



upper spectrum shows once more the part of the rotation- 
ally resolved fluorescence excitation spectrum recorded 
without the use of the selector, i.e., the free-flight spec- 
trum, shown already in [Figure 2") The spectrum in the 
middle trace demonstrates that molecules in the Oqo, Iqi 
and 2o2 rotational levels are not transmitted at all when 
the selector is operated at an ac switching frequency of 
2.70 kHz and a duty cycle of 0.50. This switching fre- 
quency is below the abrupt low-frequency cut-off shown 
in Figure 5 and it is evident from the spectrum that only 



molecules in higher lying rotational levels are transmitted 
under these conditions. To create a beam containing only 
ground-state benzonitrile molecules, the selector has to 
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FIG. 7. Effective dipole moments, peff, of low rotational quan- 
tum states of benzonitrile, labeled as JRaKc^: as a func- 
tion of the electric field strength. The dashed horizontal line 
indicates the 4.52 Debye permanent electric dipole moment 
of benzonitrile. The range of electric field strengths experi- 
enced by the molecules while in the selector is indicated by 
the shaded area. In the inset the field-free energies of the 
nine rotational levels probed in the measurements shown in 
[Figure 2] are depicted. 



be operated at ac switching frequencies above those of the 
transmission maximum, under conditions that the trans- 
mission curve has a rather steep frequency cut-off on the 
high-frequency side as well. This has been done in the 
recording of the lower spectrum, during which the selec- 
tor was operated at an ac switching frequency of 3.8 kHz 
and at a duty cycle of 0.565. In this case, only molecules 
in the Oqo and in the 2o2 rotational levels as well as a very 
small fraction of molecules in the Iqi rotational level are 
transmitted. 



IV. DISCUSSION 

To quantitatively describe the operation characteris- 
tics of the selector we have to look in detail to the inter- 
action of benzonitrile in any of the low-lying rotational 
levels with electric fields. In the following, we will re- 
strict ourselves to the Oqo, the Iqi and the 2o2 levels, 
the lowest ones of the nine rotational levels probed in 
[Figure 2\ whose field-free energies are schematically de- 
picted in the inset to [Figure 7| In general, the force 
that neutral polar molecules experience in the electric 
field is given by the negative gradient of the Stark en- 
ergy of the quantum level that they are in. This force 
can also be expressed as the product of the gradient of 
the electric field strength with an effective dipole mo- 
ment, Peff- This effective dipole moment is the negative 
derivative of the Stark energy of the particular quantum 
state with respect to the strength of the electric field. 
In the presence of an electric field, every JxaK^ level is 



split into J + 1 different M-components, of which M = 
is non-degenerate whereas all the M 7^ components 
are doubly degenerate. The calculated effective dipole 
moments of the individual M components of the three 
rotational levels mentioned above are shown as function 



of electric field strength in [Figure 7| At low electric field 
strengths there are several avoided crossings between lev- 
els of equal symmetry, resulting in a complicated depen- 
dence of Peff on the strength of the electric field. At 
high electric field strengths, the effective dipole moments 
all converge to a value of 0.0759 cm~^/(kV/cm), cor- 
responding to the permanent electric dipole moment of 
benzonitrile of 4.52 Debye. In the selector, the molecules 
that are on stable trajectories experience electric fields 
in the range indicated by the shaded area in [Figure 7| 
For benzonitrile molecules in the lowest rotational level, 
in particular, this implies that the value for the effective 
dipole moment is almost constant. It is also seen that 
the effective dipole moment is the largest for the ground- 
state level, explaining why it is completely absent from 
the transmission spectrum shown in the middle panel of 



Figure 6 



The effective dipole moment of the Jk^k^M = 2o22 
level is indistinguishable from that of the Oqo ground- 
state level in the electric fields inside the selector. The 
value of Peff foi" the M = 1 and M = components of 
the 2o2 rotational level are lower, but are in turn almost 
identical to those of the M = 1 and M = components 
of the loi rotational level, respectively. This grouping of 
effective dipole moments originates from the electric field 
dependence of the Stark shifts in the so-called pendular 
or harmonic librator limit [23l |24] ; although the group- 
ing of effective dipole moments can be readily derived 
for a diatomic molecule [21] and even for a symmetric- 
top molecule [25], the correlation of the eigenenergies in 
the strong-field limit with those of the free rotor has not 
been rigorously analyzed for an asymmetric-top molecule 
like benzonitrile. Upon exiting the electric fields of the 
selector, the population in each of the M components 
adiabatically transfers to the (degenerate) field-free ro- 
tational levels. In the LIF detection zone, therefore, the 
summed population of all M components belonging to a 
certain rotational level is measured. This can nicely ex- 
plain the degree of transmission of benzonitrile molecules 
in the 2o2 level through the selector as observed in the 
lower spectrum in Figure 6 The integrated LIF signal 



intensity of ground-state molecules in this spectrum is 
seen to be decreased to about 3/4 of its value in the free- 
flight spectrum, whereas this is about 1/3 for molecules 
in the 2o2 level. Considering that 40% of the molecules 
entering the selector in the 2o2 level will be in the M = 2 
component, the observed intensity in the spectrum is as 
expected when only the ground-state level and the M = 2 
component of the 2o2 level contribute to the LIF signal 
behind the selector. From the weak LIF signal of ben- 
zonitrile molecules in the Iqi level in this spectrum it 
can be deduced that molecules in the M = 1 component 
of this level (as well as of the 2o2 level) have at least 
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FIG. 8. Upper row: Measured intensity of the LIF signal, 
plotted in a false-colour representation, as a function of the 
velocity in the x-direction and as a function of the arrival 
time for two different ac switching frequencies. Data from 
[Figure 4] Lower row: Corresponding simulated arrival time 
distributions. The calculated velocity distributions are con- 
voluted with a Lorentzian profile with a FWHM width of 
2.2 m/s. 



a factor five lower transmission than the ground-state 
level. This nicely demonstrates the degree of selectiv- 
ity that can be obtained with the present selector as it 



is seen from Figure 7 that the relative difference in the 
"electric-field- aver aged" effective dipole moment of this 
M = 1 component and the ground-state level is less than 
10%. 



In Figure 8 the experimental data shown in Figure 4 
are compared to the outcome of trajectory calculations. 
In these calculations a pulse of molecules with an ini- 
tial length of 4.0 cm (FWHM; Gaussian distribution), 
with a mean forward velocity of 570 m/s and a veloc- 
ity spread of 25 m/s (FWHM; Gaussian distribution) is 
assumed. The initial transverse velocity distribution is 
described by a Gaussian distribution with a FWHM of 
10 m/s, centered around zero velocity. The dimensions 
of the molecular beam machine and the characteristics of 
the switching sequence as applied in the experiment are 
used as input. Effects due to the fringe fields near the 
ends of the electrodes of the selector are neglected in the 
simulations. The electric fields in the selector are calcu- 
lated using a finite element program package (COMSOL 
Multiphysics) and are assumed to be the same over the 



full 1640 mm length of the selector. Switching between 
the two electric field configurations is assumed to be in- 
stantaneous. The voltages applied to the electrodes have 
been set to ±9.5 kV, just as used in the experiment. The 
electric field dependent dipole moment for the Oqo level 
as shown in [Figure 7| is taken, and the equations of mo- 
tion are solved in sufficiently small time-steps. In the 
detection region, molecules in a 2.5 mm high interval are 
accepted, weighed by the duration that they are in the 
2.5 mm diameter cylindrical interaction volume, i.e., as- 
suming that the laser excitation step is not saturated, 
and the number of molecules are plotted as a function 
of their arrival time (binned in 2 ps intervals) and as a 
function of their velocity component in the x-direction 
(binned in 0.3 m/s intervals). To enable a direct com- 
parison with the experimentally observed patterns, the 
calculated velocity distributions at a given arrival time 
are convoluted with a Lorentzian profile with a FWHM 
width of 2.2 m/s, thus taking the 8 MHz homogeneous 
broadening due to the finite lifetime of the electronically 
excited state, that is inherent in the measurements, into 



account. It is seen from Figure 8 that the experimentally 
observed patterns are quantitatively reproduced in the 
simulations, implying that in particular also the trans- 
verse velocity distributions of the molecules exiting the 
selector are described correctly. 



The calculated time-integrated and time-resolved num- 
ber of ground-state benzonitrile molecules transmitted 
through the selector as a function of ac switching fre- 

Also in these calcula- 



quency is shown in Figure 9 



tions, the velocity distributions along the x-direction at a 
given arrival time are first convoluted with a Lorentzian 
profile with a FWHM width of 2.2 m/s and then only 
the molecules within a 0.5 m/s interval centered around 
m/s are accepted, thereby mimicking the LIF detection 
of the molecules in the actual experiment and enabling 



a direct comparison to the measurements shown in Fig- 
ure 5 1 In the comparison of the calculated and measured 



transmission curves shown in the upper row, it should 
be noted that the vertical scale of either set of curves, 
shown on the right and the left, respectively, is normal- 
ized to the situation in free flight. The observed overall 
transmission of the selector is seen to be about a factor 
of three lower than expected from the calculations. The 
shapes of the measured and the calculated curves are 
in good agreement, however. The observed modulations 
in the transmission curves and the low-frequency onset 
of the transmission of the selector are well reproduced 
in the calculations. At the high frequency end the ex- 
perimentally observed transmission curves drop slightly 
faster than expected from the simulations. The observed 
change in the transmission characteristics with changing 
duty cycle is in quantitative agreement with the outcome 
of the trajectory calculations. The latter can therefore 
be used to quantify the p/Ap resolution that can be ob- 
tained with the present selector. 
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FIG. 9. Upper row: Calculated time-integrated transmission 
curves for three different duty cycles (solid curves), compared 
to the corresponding measured curves (dashed curves), al- 
ready shown in [Figure 7] The vertical scale for the measured 
curves is on the left, while the one for the calculated ones is 
on the right. Lower row: Calculated number of ground-state 
benzonitrile molecules transmitted through the selector, plot- 
ted in a false-colour representation, as a function of the ac 
switching frequency and as a function of arrival time for three 
different duty cycles. The panel on the left shows the time- 
integrated signals as a function of ac switching frequency, also 
shown in the upper row. See text for further details. 



V. CONCLUSIONS 

In this work, we have demonstrated and character- 
ized the performance of a state-of-the-art m/p selector 
for large molecules. Using a pulsed beam of benzonitrile 
molecules in combination with detection with a narrow- 
band continuous wave laser, the quantum-state-specific 
transmission of the selector as well as the transverse ve- 
locity distribution of the molecules exiting the selector 
have been measured as a function of ac switching fre- 
quency. It has been demonstrated that by controlling 
the applied ac waveforms, a p/Ap resolution can be ob- 
tained that is sufficient to exclusively transmit molecules 
in quantum states that have the same m/p- value as the 
absolute ground state of benzonitrile. The measured 



transmission characteristics of the m/p selector are in 
quantitative agreement with the outcome of trajectory 
simulations; only at high frequencies the observed trans- 
mission is smaller than expected. It is well-known, how- 
ever, that the transmission through the selector under 
these operating conditions is extremely sensitive to mis- 
alignments [13l|2T]; fluctuations in the voltages applied to 
the electrodes might also play a role. Moreover, at these 
high frequencies the minute losses due to nonadiabatic 
transitions that might occur each time that the fields are 
being switched will be most notable. If we assume, for 
instance, that some 5% of the ground-state molecules are 
lost to other quantum states each time that the fields are 
being switched while the molecules are inside the selec- 
tor, we can bring both the observed overall transmission 
efficiency and the observed drop in transmission at high 
ac switching frequencies in perfect agreement with the 
calculations. 

The good agreement between the observed and cal- 
culated performance of the m/p selector is not only a 
testimony to the accurate mechanical design of the selec- 
tor, but it also implies that trajectory calculations can 
be used to predict the outcome of future experiments. In 
the present selector, the number of full switching cycles 
that the molecules experience near optimum transmis- 
sion is only about ten, which poses an intrinsic limita- 
tion to the attainable p/Ap resolution. Nevertheless, tra- 
jectory calculations predict that for molecules with the 
mass of benzonitrile and with an effective dipole moment 
of around 4.5 Debye, the selector can be operated such 
that the overall transmission for levels with a constant 
effective dipole moment differing by 5% can be made to 
differ by a factor of two; when the p/Ap resolution is 
defined in this way, a value of 20 is obtained. If a beam 
containing mainly the species with the lower dipole mo- 
ment of the two is desired, the selector can operate be- 
low the low-frequency cut-off for the other species at a 
duty cycle of 0.50. A beam predominantly containing 
the species with the largest dipole moment of the two 
can only be obtained at higher frequencies and at higher 
(or lower) duty cycles, thereby necessarily reducing the 
overall transmission somewhat. Apart from separating 
quantum states of a certain molecule, the selector can 
also be used to separate molecules with different masses. 
It is expected, for instance, that in the current setup 
benzonitrile-Ar clusters will be present in the molecu- 
lar beam as well. These clusters are expected to have 
a dipole moment that is very similar to that of the bare 
benzonitrile molecule, but their rotational level structure 
will be more dense and their mass is considerably higher. 
With the selector operating at a frequency below the low- 
frequency cut-off for benzonitrile, with a duty cycle of 
0.50, it should be possible, therefore, to selectively trans- 
mit these higher mass clusters. A further possibility to 
increase the p/Ap resolution might be the modulation of 
the high voltages applied to the electrodes with specific 
frequencies to selectively eject molecules with unwanted 
m/p ratios via parametric amplification of their trans- 
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verse motion. Studies are underway to investigate these 
applications of the m/p selector. 

The m/p selector provides samples of highly oriented 
molecules. To maintain a defined orientation for a subse- 
quent experiment, a static electric guiding field of about 
25 kV/cm is needed at the exit of the selector, which is 
well feasible. It should be mentioned that as an alterna- 
tive to the m/p selector, a static "two- wire" deflector has 
been used to prepare pure samples of large molecules for 
orientation and alignment studies [26, 27 as well as to 
separate structural isomers of molecules [28]. Although 
the implementation of the deflector is more straightfor- 



ward than that of the m/p selector, its application is also 
more limited; the difference in applicability between the 
two devices is similar to the difference in applicability be- 
tween a static ion deflector (a separator or bender) and 
a quadrupole mass- filter. 

VI. ACKNOWLEDGMENTS 

The design of the electronics by G. Heyne, V. 
Platschkowski and P. Schlecht has been crucial for this 
work. We acknowledge the software support by U. Hoppe 
and H. Junkes. 



[1] G. Scoles, ed., Atomic and molecular beam methods, Vol. 
1 & 2 (Oxford University Press, New York, NY, USA, 

1988 & 1992) ISBN 0195042808.^ 

D. H. Levy, Science, 214, 263 (198l 7 
T. R. Rizzo, Y. D. Park, L. Pe teanu, and D. H. Levy,|X| 



Chem. Phys., 83, 4819 (1985) 



S. Y. T. van de Meerakker, H. L. Bethlem, and G. Meijer, 
Nature Phys., 4, 595 (2008) 

J. P. Gordon, H. J. Zeiger, and C. H. Townes, |Phys.| 
Rev., 99, 1264 (1955) 

D. Auerbach, E. E. A. Bromb erg, and L. Wharton, [X] 
Chem. Phys., 45, 2160 ( 1966)] ^ 

E. D. Courant and H. S. Snyder, Ann. Phys., 3, 1 (1958), 
D. Kakati and D. C. Laine, Phys. Lett. A, 24, 676 (1967) 
H. L. Bethlem, A. J. A. van Roij, R. T. Jongma, and 
G. Meijer, Ph ys. Rev. Lett., 88, 133003 (2002) 

J. van Veldho ven, H. L. Bethl em, and G. Meijer, |Phys.| 



Rev. Lett., 94, 083001 (2005) 



T. Junglen, T. Rieger, S. A. Rangwala, P. W. H. Pinkse, 
and G. Rempe, Phys. Rev. Lett., 92, 223001 (2004) 
T. E. Wall, S. Armitage, J. J. Hudson, B. E. Sauer, J. M. 
Dyne, E. A. Hinds, and M. R. Tarbutt, |Phys. Rev. A,| 
"80, 043407 (2009), 



K. Wohlfart, F. Gratz, F. Filsinger, H. Haak, G. Mei jer, 
and J. Kiipper, Phys. Rev. A, 77, 031404(R) (2008)' 

F. Filsinger, U. Erlekam, G. von Helden, J. Kiipper, and 

G. Meijer, Phys. Rev. Lett., 100, 133003 (2008), 

F. Filsinger, S. Putzke, H. Haak, G. Meijer, and 
J. Kiipper, Ph ys. Rev. A, 8 2, 052513 (2010) 
K. Wohlfart, M. Schnell, J.-U. Grabow, and J. Kiipper, 
|J. Mol. Spec, 247, 119 (2008)| 



[is: 

[19 

[20; 

[21 



[22; 

[23 
[24 

r 

[25; 

126 



[27 
[28 



D. R. Borst, T. M. Korter, and D. W. Pratt, |Chem. 
Phys. Lett., 350, 485 (2001) 

C. M. Western, "PGOPHER, a program for simulating 
rotational structure," University of Bristol, Bristol, UK, 
URL: http : / / pgopher . chm . br is . ac .ukj 
T. Kobayashi, K. Honma, O. Kajimoto, and S. Tsuchiya, 
J. Chem. Phys., 86, 1111 (1987) 

F. Lahmani, C. Lardeux-Dedonder, and A. Zehnacker, 
Laser Chemistry 10, 41 (1989) 

H. L. Bethlem, M. R. Tarbutt, J. Kiippe r, D. Carty, 
K. Wo hlfart, E. A. Hinds, and G. Meijer, |J. Phys. B,| 
39, R2 63 (2006) 

S. Schlunk, A. Marian, P. Ceng, A. P. Mosk, G. Mei- 



jer, and W. Schollkopf, Phys. Rev. Lett., 98 (2 007)] [doT 
10.1103/PhysRevLett.98.223002 



_H. J. Loesch and A. Remscheid, J. Chem. Phys., 93, 4779 
(1990) 



B. Friedrich and D. R. Herschbach, |Z. Phys. D, 18, 153 
(1991) 



M. Hartelt and B. Friedrich, The Journal of Chemical 
Physics, 128, 224313 (2008) 
L. Holmegaard, J. H. Nielsen, I. Nevo, H. Stapelfeldt, 
F. Filsinger, J. Kiipp er, and G. Meijer, Phys. Rev. Lett.,| 



102, 023001 (2009) 



F. Filsinger, J. Kiipper, G. Meijer, L. Holmegaard, J. H. 
Nielsen, I. Nevo, J. L. Hansen, and H. Stapelfeldt, pT| 
Chem. Phys., 131, 064309 (2009) 

F. Filsinger, J. Kiipper, G. Meijer, J. L. Hansen, J. Mau- 
rer, J. H. Nielsen, L. Holmegaard, and H. Stapelfeldt, 
|Angew. Chenr7l21, 7033 (2009), 



